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Abstract 

^ ■ We study the analytical continuation in the complex plane of free 

•/^ . energy of the Ising model on diamond-like hierarchical lattices. It is 

I known [|l2| , 13 1 that the singularities of free energy of this model lie on 

the Julia set of some rational endomorphism / related to the action 
, of the Migdal- Kadanoff renorm-group. We study the asymptotics of 

free energy when temperature goes along hyperbolic geodesies to the 
boundary of an attractive basin of /. We prove that for almost all 
(with respect to the harmonic measure) geodesies the complex critical 
exponent is common, and compute it. 
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1 Introduction 

The purpose of this article is to analyse complex singularities in temperature 
of the free energy J-' in the Ising model on diamond-like hierarchical lattices. 
According to the traditional point of view a phase transition manifests itself 
as a singularity of as a function of thermodynamic parameters (like tem- 
perature and external magnetic field). From this point of view the theory 
of phase transitions should describe the domain of analyticity of JF and the 
type of its singularities at points of phase transition (see [1], where diverse 
approaches to the first of these problems are discussed). 
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Since T is real analytic outside of points of phase transition, it can be 
continued into complex space with respect to the thermodynamic param- 
eters. Description of its complex singularities is of great interest for the 
theory of phase transitions because it determines analytic properties of the 
thermodynamic function . 

The celebrated Lcc-Yang theory (sec [2]) gives a realisation of this ap- 
proach describing the singularities of the analytic continuation of the free 
energy in the ferromagnetic Ising model with respect to the external mag- 
netic field. It proves that the zeroes of the grand partition function in the 
ferromagnetic Ising model lie on the imaginary axis, and hence complex sin- 
gularities of the free energy lie on the imaginary axis as well. An important 
problem stated in [2] is to study the limit distribution of zeros of the grand 
partition function, since the free energy can be expressed as a logarithmic 
potential over this distribution. 

The problem of description of complex singularities of the analytic con- 
tinuation of thermodynamic functions in temperature is also very interesting 
from different points of view. Many properties of asymptotic behavior of ther- 
modynamic functions in vicinity of a critical point were investigated through 
the Kadanoff-Wilson-Fisher renormalisation group theory (see e.g. [3,4]). 
It gives a local form of real critical singularities of thermodynamic functions 
which has a nice universal scaling structure. A problem is how are these 
singularities continued to complex space and what is their global structure 
in complex space? 

Unfortunately no general theory like the Lee-Yang theory exists which 
describes for general models global complex singularities of thermodynamic 
functions in complex temperature plane . However some exact results were 
obtained for the two dimensional Ising model. The main tool here is the fa- 
mous Onsager solution. It turns out that in isotropic two dimensional Ising 
model the zeroes of the partition function lie asymptotically on two circles 

g-2J/T = ±1 + ^giv' (it 

was conjectured by Fisher [5] and proved in [6]). 
Later it was shown (see [7,8] and references there) that in anisotropic two 
dimensional Ising models on diverse lattices the zeroes of the partition func- 
tion fill some planar regions in the complex temperature plane and in some 
cases the density of the limit distribution of zeroes can be found explicitly. 

In the present paper we consider another exactly solvable model, namely 
the Ising model on diamond-like hierarchical lattices (see [9-11]) . In this 
case the Migdal-Kadanoff renormalization transformation turns out to be to 
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a rational map / on the Riemann sphere which can be found explicitly. The 
following nice observation was made in papers [12, 13]: the set of complex 
points of phase transition coincides with the Julia set of this map. In this 
paper we study the analytical properties of the free energy F near the sin- 
gular points. We believe that clear understanding of analytical properties 
of thcrmodynamical functions in the complex plane has much to do with 
physical nature of the model. 

The hierarchical sequence of the diamond-like lattices depends on one 
natural parameter b > 2. The lattice To is just two "outer" sites related by 
a bond. In order to obtain Fi we insert in between the outer sites b inner 
sites related by bonds with outer ones (see Fig. for 6 = 3). Then in order 
to obtain F„+i we replace each bond in F„ by the lattice Fi, n = 1, 2, .... 




Figure 



We will refer to papers [11-15] for the description of the Ising model on 
these lattices and the calculation of the thcrmodynamical functions. The 
starting point for us is the following explicit formula for the free energy: 

T rp OO 

•^=-9 - 9 E(2&)-"Ml + <) (1-1) 

^ ^ n=0 
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where b = 2*^"^, d is the "dimension" of the lattice, J is the interaction 
constant, T is temperature, and the sequence tn, n > 0, is given by the 
following recurrent equation 

tn+l^f{tn), n>0 (1.2) 

where 

m = (1.3) 

with the initial data 

to = exp(-^)^G(T). (1.4) 

Equations (1.2)- (1.4) mean from the physical point of view that the map 
T G^^ o / o G{T) gives the rescaling of temperature under the Migdal- 
Kadanoff renorm-group transformation (see [9-15]). Note that the points 
t — and t — 1 are superstable fixed points of the map t ^ f{t) (low- and 
high-temperature fixed points of the renorm-group) and for fo > 1 there exists 
the unique unstable fixed point tc on [0,1] . The critical temperature T^, is 
equal to G~^(tc)- 

Formulas (1.1)-(1.4) make sense for complex values of T as well. So, we 
can consider the analytical continuation of free energy JF(T) from the positive 
axis T > into the complex plane. It is not hard to see that the singularities 
of JF lie on the Juha set J(/) [12, 13] (For the definition of the latter, see 
one of the surveys [7-10]). 

Let us consider now the immediate attracting basins Qq and Qi of points 
and 1. One can show that C \ J(/) is the union of preimages of these 
domains, and is a Jordan domain in J(/) (see §2). In this paper we study 
the boundary properties of !F in the domain Qq- To this end let us consider 
the Riemann map ip : Qq — > U of Qq onto the unit disk. The hyperbolic 
geodesies in flo are just the images of the radii in U. Denote by Bt the 
geodesic ending at r e d^o- Let us consider also the harmonic measure /i on 
dQo, i-S- 1^ — "fp*^"^: where A is the Lebesgue measure on the circle d\J=T. 
For t e Bt denote by l{t) the length of B^- from t to r (perhaps, l{t) — oo). 

In the present paper we prove: 

(i) The derivative J-'' of free energy is continuous up to the boundary of 
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(ii) For b > 2 the second derivative is discontinuous in Q, and has the 
following asymptotics on ^-almost all geodesies: 



In t) ln2 ^ 1 ,^ , 

hm — ' , ,; i ' =ac^l- — ^l- - — -. (1.5) 

t^T,tGB^ -\nl{t) In 6 d-1 ^ ^ 

This means that for almost all geodesies the specific heat critical exponent 
in the region of low temperatures is universal and equal to 1 — 

Now let us dwell in more detail on the content of the paper. 

In §2 we describe the dynamical properties of the endomorphism /. In 
particular, we show using the Douady-Hubbard straightening theorem that 
Qq is a Jordan domain. 

In §3 we show that J-'' is continuous in clflo and that dflg is the natural 
boundary of analyticity of J^. The proof is based upon some amusing obser- 
vations concerning / (its relation to the Koebe function and a Tchebyshev 
polynomial) . 

In §4 we discuss some technical background: the Bowen-Ruelle-Sinai ther- 
modynamical formalism and the construction of the natural extension (the 
inverse hmit) of /. These are the main tools (together with the ergodic 
theorem) for the accurate computation of the critical exponent. 

In §5 we discuss the functional equation for JF" and related spectral prop- 
erties of the weighted substitution operator in the disk-algebra. 

§6 is the central section of the paper: here we give the computation of 
the critical exponent, provided T" is not continuous up to the boundary of 

In §7 we prove that J-'" really satisfies this property, which completes the 
proof of the main result. 

In the last §8 we discuss some related problems. 

Acknowledgement. We are grateful to J. Milnor for looking through 
the text and making useful remarks, to M. Fisher and the referee for critical 
comments yielding the improvement of the exposition, to J. Milnor and G. 
Tusnagy for making nice computer pictures, and to NEFIM fund of Hungar- 
ian Academy of Sciences for the support of the visit of one of the authors 
(P.M.B.) to Budapest. 
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2 Dynamics of the map / : t 



We refer to the surveys [16-19] for the general view of the dynamics of complex 
rational maps. We will use some concepts and facts of this theory without 
extra explanations. 

Let us introduce the following notations: 

yon _ J Q Q J jg ^Yic n-fold iterate of /; 

C{f) is the set of its critical points (a rational map of degree d has 2d — 2 
critical points counting with multiplicity); 
1 1 ■ 1 1 is the spherical metric on C; 
U — {z :\ z \< 1} is the closed disk; 
U° = intU is its interior; 
T = (9U is the unit circle; 
B{a,r) = {2; e C :| ^ - a |< i?} for a e C; 
J(/) is the Juha set of /. 
The function 



is related to the well-known extremal Koebe function{see [20]) 

}Co{z) = - 



Setting K{z) = -4/Co(-^) and S{z) = Sb{z) = z^ we have f^KoS. 

The relation of / to the Koebe function is quite mysterious, especially if 
one relates the coefficient 4 to the Koebe constant 1/4. It becomes still more 
amusing if to observe that K{t) is conformally conjugated to the Tchebyshev 
polynomial T : t ^ 2t^ — 1. Indeed, the function K has two simple critical 
points Cx — \ and C2 = —1. Moreover, /(ci) = Ci, i.e. C\ is superstable fixed 
point, and C2 i-^ cxd t-^ 0, where is repelling fixed point. Up to conformal 
conjugation, T is the unique rational function of degree 2 possessing such 
properties. More specifically, ipoKo{p~^ = T where ip : t ^ is the Mobius 
transformation mapping the triple {0, 1, 00} onto the triple {1, cxd, — 1}. In 
particular, it follows that the Juha set J{K) coincides with the negative 
semi-axis [—00, 0] = (f~^[—l, 1]. 

The power functions and Tchebyshev polynomials play a particular role 
in the iteration theory. They appear as the exceptions in a number of prob- 
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lems; e.g., only these functions have a Juha set with simple geometry. The 
composition f = K o S does not possess such a property (see Fig. 1). 

A rational function g is called critically finite if the orbits {g"'{ci)} of all 
its critical points are finite. 

By the chain rule 

Cif) = CiS) U S-'CiK) = {0, 00} U {a,}ti U {/3,}5=i, 

where are the 6th roots of 1 , and f3j are the 6th roots of -1. Moreover, 
and ai = 1 are superstable fixed points which absorb the orbits of all other 
critical points: 1, /5j t-^ 00 1-^ 0. Thus, the function / is critically finite. 

Denote by Qa = ^{a) , the component of F{f) containing a. The domains 
Qq and Qi are called the immediate basins of the fixed points and 1. 

We say that a rational function g satisfies the axiom A if the following 
equivalent properties hold : 

(i) The orbits of all critical points converge to stable cycles; 

(ii) g is expanding on the Julia set, i.e., there exist constants C > and 
A > 1 such that 

II dg'^iz) \\> CA"(^ G J(/),n G N). 

It follows from above that our function / satisfies (i) and hence satis- 
fies axiom A. This implies in particular that the Fatou set consists of the 
preimages of the immediate basins Qq and Qi. 

Set Q = c\Qo and T = dflo (these notations will be used up to the end of 
the paper). We will show now that F is a Jordan curve and even a quasicircleQ 
(but by Fatou's theorem it has no tangents at any point). To this end we 
apply the Douady-Hubbard straightening theorem (see [21]). 

Let V and V be two simply connected domains bounded by piecewise- 
smooth curves, and cl\^ G V G C. A map g : V ^ V is called polynomial- 
like of degree d if it is a ci-sheeted analytical covering of V over V having 
no critical points on dV. By the Riemann-Hurwitz formula, such a map has 
d — 1 critical points in V counting with multiplicity. Set 

K{g) = {z : g^{z) G V{n = 0, 1, ...)}, K\g) = mtK{g). 

The K{g) is a compact subset of V. 
-'^ actually, the last holds automatically 
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The straightening theorem states that any polynomial-hke map g is quasi- 
conformally conjugated to a polynomial of the same degree, i.e., there exists a 
quasi-conformal homeomorphism ^0 : C — > C such that ipog \ W — hoip \ W 
for some domain W, K gW gV. Moreover, ip \ K'^{g) is conformal and 

il}{K{g)) = C\{z : h°'\z) ^ oo (n ^ oo)} 

is the filled-in Julia set of h. 

Lemma 2.1 The domain is Jordan, and its boundary is a quasi-circle. 
The restriction / | JIq is conformally conjugated to the power transformation 
of the unit disk U°. 

Proof. Let us construct a neighbourhood of Q on which / is a polynomial- 
hke map of degree b. To this end note that the function K conformally maps 
the disk U'^ onto the plane slitted along the semi- axis Ri = [1, oo) (this is 
the characteristic property of the 4-fold Koebe function). Let us consider the 
domain V{see Fig. 2) bounded by the arc 71 of the circle B{l,e), the arc 72 
of the circle B{0,K) where i? >> 1 and two horizontal intervals. 




Figure 2 
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Let V be the component of the inverse image / ^{V) containing 0. Then 
clV C U° since f{dlJ) = Ri lies outside V. Besides, dV n 71 = for 
sufficiently small e. Indeed, as 1 is a stable fixed point, the arc /(71) hes 
inside the disk B{l,e) and, hence, outside V. 

Thus, clV C v. 

Further, it is clear from V = {K \ U)^^ o [S^^V) that V is simply- 
connected. Indeed, it is elementary that S^^V is simply-connected (see e.g., 
[17], Lemma 1.4), while X | U is univalent. 

We have shown that f : V ^ V is a, polynomial-like map. Its degree 
is equal to h since V <Z Uq contains the unique [h — l)-fold critical point 0. 
Clearly ^ C K{f \ V). 

By the Straightening Theorem, / : ^ is quasi-conformally conju- 
gated to a polynomial h of degree h. Normalize the conjugating homeomor- 
phism if) in such a way that '?/'(0) = 0- Then becomes a (6 — l)-fold critical 
point for h. It follows that h{z) = Cz^. Normalizing ip additionally in such 
a way that ip{tc) = 1 where tc is a real fixed point lying on F ,we get C = 1. 

Thus, ip conjugates f : V —>■ V to the power polynomial h : z ^ z''. 
Consequently, = '0~^(U) is a Jordan domain bounded by a quasi-circle, 
and ip conformally conjugates / | JIq to /i | U. The lemma is proved. g 

Remcirk. Tan Lei showed us another proof of the above lemma which 
can be applied to the high-temperature region as well. 



3 Analytic properties of J^' 

In this section we will show that the derivative J^'(t) is continuous in the 
closed set fl. So, J-'' belongs to the disk-algebra A{fl), i.e. the algebra of 
functions continuous in D, and holomorphic in JIq- We will get this estimating 
the derivative || Df W^, in a special Riemann metric /i. 

From now on we will consider the following function F instead of free 
energy (1-1): 

where g{t) — ln(l +t^). Clearly, its analytical properties are the same as 
those of J^. Note that g is analytic in a neighbourhood of Q, and so series 
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(3.1) converges uniformly in Q. Hence F e A{Q). Further, for z & Qq we 
have 

CO 1 

n^) = E7^^'(r(^))(/7(^) (3.2) 

n=o y"^^) 

We want to show that this series converges uniformly in Q, which cer- 
tainly implies F' e A{fl). The required statement follows from the following 
estimate: 

Lemma 3.1 | (/")'(^) |< C(V26)" for z e n. 

Proof Let us recall that f = K o S (see §2). The power function 5" satisfies 
the functional equation exp(te) = ^(cxp^;). From the dynamical viewpoint 
it means that exp semiconjugates the transformations L : z ^ hz and S. 
Denote by a the Euclidean metric on C, and by = exp^ a its image on the 
punctured plane C* = C\{0}. We have \ diJL\ = \ dt\ / \ t |. 
As II DL{z) ||o.= 6 for 2; e C, 



DS{t)\\^^h (3.3) 



for i e C*. Besides 



where i^{t) — {1 -\- t) / {1 — t) . It is surprising that exactly this function 

conjugates K to the Tchcbyshcv polynomial T : t ^ 2r^ — 1. Due to this 
observation it is reasonable to pass to the conjugated function 

where it! = (f) o S o (f)"^. Consider the corresponding Riemannian metric 
= By (3.3), (3.4), 



Hence for r = (pit) 



DR{T)\l^b, ||L'r(T)||,= 



Df{t)\\,^\\Dh{r)\l=y^^^. (3.5) 
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The function h has a superstable fixed point 1 = (/^(O) with the immediate 
attracting basin = ipQo. Since Q C U°, the set W = clW^ = ip{Q) lies 
in the right half-plane (p{\J^) — — {r : Rer > 0}. As W is g'-invariant, 
T{R{W)) C P° and hence (see Fig. 3) 

R{W) C T-^(P°) = {r = x + iy : -y^ > 1/V2} 
Consequently, | R{t) \> l/\/2 for r e W, and (3.5) implies 

\\Df(t)\\^<bV2, ten 

Hence, 

\\Dr{t)\\^<{bV2r, neN. 

But on the boundary dn the metric is equivalent to the Euclidean metric 
and, hence, | (/")'(t) |< (7(^26)" for t e dn. By the Maximum Principle, 
this inequality holds for ten. The lemma is proved. g 




Figure 3 



Now let us establish some global analytical properties of F'. It is curious 
that this function has no singularities at points Pi = \[—\. 

Lemma 3.2 (i) The junction F' is a single-valued holomorphic function on 

the Fatou set N{f). 

{u)The set 0° is the maximal domain of analyticity of F'. 
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Proof. Consider the multy- valued function cr{t) = g{t) + -^g{f{t)) where 
g{t) = ln(l+t'') as above. Let us show that it is regular near Pj = j = 

1, b. We have 

.W=ln(l+t') + ^ln(l + ^^) = 



2b 

So, cr(/5j) = iln4 + (I + 7in)i, and we see that Pj are regular points for 
a. Hence, a is regular in the components fi{Pj)- 

Now let V be an n-fold preimage of some ^{Pj)- By (3.1), 

F{t) — (^2b)^ ^^^^^^ ~'~ regular function. 

Consequently, F is regular in V. 

Thus, n(oo) is the only component in which F is not regular, and there 
we have 

F{t) = ln(l + t'') + regular function. 

Hence, 

F'{t) = j^^i, + regular function 

is regular in ^(oo). 

So, F' is regular on the whole Fatou set A/"(/),and it is obvious from (3.2) 
that it is single- valued. The (i) is proved. 

(ii) Let us consider the functional equation for F: 

l_Fof-F = g (3.6) 
Taking the derivative, we get 



Lp'of-F'^g'^^—^. (3.7) 



Provided F' can be analytically continued beyond Vf' into some neighbour- 
hood t/ of t e dVL , it follows from (3.7) that it can be continued as a 
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meromorphic function into /"f/. Since /"f/ D J{f) for some n , the function 
F' is meromorphic on the whole sphere, i.e. rationaL 

But we have shown that F' has no poles in N{f). If F' had a pole at 
ate J(/) then by (3.7) it would have had poles at all points of the grand 
orbit 

oo oo 

o{t) = u U mrt) 

m=0 n=0 

But it is impossible since this orbit is infinite. So, F' has no poles at all, and 
this absurdity completes the proof. a 



4 Symbolic dynamics, thermodynamical for- 
malism and the natural extension of / | Q 

To study the boundary behaviour of F" we need the Bowen-Ruelle-Sinai 
thermodynamical formalism (see[22],[24] or [18]). Wc will state the main 
results of this theory for our particular map / | F. The theory can be 
applied to this map without any problem because it is expanding (sec §2). 

Let us consider the homeomorphism : Q ^ \J conformal on Qq and 
reducing f \ Q to S : z i-^ z'^ (Lemma 2.1). Let ip{t) = re^"^*^ e U. Associate 
to a point t & ^ the pair (r, F+) where r G [0, 1] , and = {so,Si, ...) is the 
6-adic decomposition of the 6 e [0, 1). Then / turns into the transformation 
(r, £+) 1-^ (r'', cr_|_£+) where 

(J+ : (£o,£i, ■■■) ^ (£i,£2, ■■■) 

is the shift on the space of all one-sided 6-adic sequences. Sometimes we 
will identify Q with [0, 1] x T,^ , though the described correspondence is not 
one-to-one. Then / | F will be identified with the shift 

Denote by [eo,...en-i] 6-adic cylinders in and corresponding 6-adic 
intervals in F. Let Bg — ^~^{re^'^*^ : < r < 1} be hyperbolic geodesies in 
f2, and F^ = ^"^{re^'^'^ : < ^ < 1} be equipotential levels (for the Green 
function In | ip |). 

Let p be a Holder function on F which is called the potential (here we 
pass from electrostatics to thermodynamics). Set 

k=0 
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The Gibbs measure Up on F corresponding to the potential p is the measure 
satisfying the following estimates on cylinders: 



PpleQ-.-Sn-i] X exp[S'„p(tso,„s„_J - nP], (4.1) 

where is any point of [£:o---£^n-i], and P = Pf{p) is a constant called 

the pressure, and the sign "ax /5" means CiP < a < C2/5. 

The main result of the Bowen-Ruelle-Sinai theory states that for any 
Holder function p there exists the unique Gibbs measure z/p. This measure 
satisfies the Variational Principle 

sup {h^if) + I fdu) = K^{f) + / fdup = Pfip) (4.2) 

where M{f) denotes the compactum of all /-invariant probability measures 
on r, and h^{f) is the entropy of u. 

The pressure Pf{p) is the smooth convex functional of p, and its differ- 
ential at p is the Gibbs measure Up (see [23]): 

U=o= J adi^p (4.3) 

To the potential p = corresponds the unique measure of maximal en- 
tropy /i = z/q, namely the BernouUy measure with b equal states (in view 
of the model S^). The entropy of this measure is equal to the topological 
entropy of / | T : /i^(/) = h{f) = In 6. 

The Riemann map ip : Q ^ U transforms p into the measure of maximal 
entropy for z ^— z^, i.e. to the Lebesgue measure on T. Hence, p coincides 
with the harmonic measure on F corresponding to (see [20]). Consequently, 

H{0) = J Hdu 

for any function H harmonic in Qq and continuous up to the boundary. 
Applying this formula to the function 

fit) I 



H{t) = In 



■^taking in account that ijj{t) ~ t, since there is the foUowing formula for the Green 
function (cf. [32]): 

In I V h lim /°" I 
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we find the characteristic exponent of /i 

X^ = J^ln\f\dfi = \nb. (4.4) 

Let us pass now to tlie crucial construction of the natural extension (or the 
inverse limit) f : Q ^ Q (see [25]). By definition, a point z G is tlie inverse 
orbit z = {zo,z_i,...), i.e. fz_^i+^_= z_i, and J : z ^ {fzo,zo,zi,...). The 
transformation / is invertible on Q, and tliere exists the natural projection 
71 : Q ^ Q, tt{z) = zo, semi conjugating / and /. All fibers of vr except zero 
one are Cantor sets. Each /-invariant measure u onT can be uniquely lifted 
to the /-invariant measure z7 on F. 

The symbolic dynamics for / generates the symbolic dynamics for /. 
Namely, Q can be identified mod with [0, 1] xEf, where = {(...e_i, eo, 
is the space of two-sided sequences. Then 

7 : {r,e) ^ {r\ae), 

where a : S;, ^ is the left shift, and Jt turns into the Bernoulli measure 
on Eft with equal states. 

The lifts Tr = Tc'^Tr of the equipotential levels will be called the solenoids, 
r = Fi. The reason is that these sets can be supplied with the structure of 
the solenoidal group T ,i.e. the inverse limit of the group endomorphism 
z z'^ of T (see [26]). Then / | F turns into a group endomorphism, and JI 
into the Haar measure on F. 

The space Q can be regarded as a continuum-sheeted Riemann surface 
over Q, the "bunch of sheets" gluing together at zero. If one cuts Q along 
the geodesic Bq, Q is folliated into the sheets L(e-) coded by one-sided 
sequences e_ = {...e_2,£i) £ ■ 

The gluing of the sheets is fulfilled by the 6-adic shift A = Ah : + 1 

where 1=(... 0,0,1), and addition is understood in the sense of the group of 
6-adic numbers . Gluing together countably many sheets corresponding to 
an orbit {yl"(£_)}J^_^ of the 6-adic shift, we get the logarithmic Riemann 
surface W{e^) f\ All inverse functions f~^{z) become single- valued on this 
surface. 

■^The inner topology of this surface differs from the topology induced from the space O 
in which it is densely immersed 
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From the dynamical point of view, Riemann surfaces W{e-) are global 
unstable manifolds of / : if t, r G W{e-), then / ^-orbits of t and r are 
exponentially drawing together. 

For a function p{t) in Q we will write p{t) = p{t), where t = 7r(t); in 
particular, | t | = | t |. Set 

-n+l 

Wsie^) = {te Wie^) :| t \> S}, S^np{t) = - ^ p{f\t)). (4.5) 

fc=0 

The drawing together of inverse orbits originating in Ws{e-) is uniformly 
exponential. It follows in a standard way that for any function Holder 
on Vt* = f2\{0}, for (5 > and t, r G WsieJ) the following estimates hold 

I 5_„p(t) - 5_„p(r) |<C(5), (4.6) 

where C{5) does not depend on n and 

5 The weighted substitution operator 

Let us consider a function (3 G AiVL) having no zeroes on F = dVt, and 
h G Aiyt). In this and the next sections we will consider the following 
functional equation: 

mUUt) - fit) = -h{t) (5.1) 

Let us consider the multiplicative cocycle Pn(t) = P{t)P{ft)...p{f"-~H) 
associated with (3. We will assume in what follows that 

(i) the cycle (3 has the positive Liapunov exponent: 

X,{f3) = J\n\(3\dp>0 (5.2) 

(ii) 3 The function In | /5 | on F is not homologous to a constant. This 
means that there are no continuous solutions ip G C(F) of the equation 

\n\(3\=(pof-f + c 
^this assumption is convenient but it is not essential 
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with any constant c. 

Motivation. By differentiating equation (3.6), we get the equation (5.1) 
for U = F" with: 

Pit) = h{t) = ^^F'{ft)f{t) + g"{t). (5.3) 

The function /3 is holomorphic in a neighbourhood of Q and has there the 
unique root at t = 0, and the function h e A{0,) (since F' e A{Q)). As- 
sumption (i) is vahd due to (4.4): 

XM(/3) = 2XM-In26 = ln^>0. (5.4) 

Assumption (ii) holds since In | /' | is not homologous to a constant on T 
[30]. B 

The equation (5.1) has the unique solution holomorphic in ■ 

oo 

U{t) = Y.Pr^m{rt), (5.5) 
n=0 

where Pn{t) is the multiplicative cocycle generated by the function /3. 

In this section we will explain why this solution is not, as a rule, contin- 
uous up to the boundary. The words "as a rule" means: for any h G A{Q) 
outside a set of first category. In §7 we will show that the concrete function 
h given by (5.3) is not excluded (so, F" is discontinuous). 

Let us consider the weighted shift operator in ^4(^2) : 

{LpU){t)^P{t)U{ft). 

Then we can rewrite (5.1) in the following way: 

{Lp - I)U = -h. (5.6) 

This leads us to the problem of spectral properties of the weighted shift 
operator. 

It is known [27] that the spectral radius of Lp in the disk-algebra can 
be calculated by the formula: 

lnr^= sup xM: (5-7) 

u&M{f) 
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and the spectrum of Lp is the unit disk: spec(L^) = {A :| A |< r^}. 

If U is an eigenfunction of Lp then the function In | /5 | is homologous to 
a constant: 

In I /3 1= In I U{ft) I - In I U{t) | +ln | A | 

which is not the case by Assumption (ii) above. So, the operator — \I is 
injective for all A. By the Banach theorem on the inverse operator (see [28]), 
for A G spec(L^) the image Im(L^ — A/) is the set of first Baire category. 
Thus, the equation LpU — XU = /i for | A |< r/3 and generic h G A{fl) has no 
solutions in A{il). 

By (5.2), rjs > 1, and hence 1 G specL^. So, the equation (5.6) is non- 
solvable in A{Q) for generic h G A{Q). For such an h, the analytical solution 
(5.4) is not continuous up to the boundary, as we have asserted. 

In conclusion let us mention a rough obstruction for (5.1) to be solvable 
related to non-invertibility of /. Let t = (to,^-i, •••) G be an inverse orbit 
of /. Iterating (5.1) we get 

n 

U{t) = P^nmit-n) - E P-k{t)h{t^k), (5.8) 

fe=l 

where P-k{t) = {P{'t-i)---P{t-k)]~^- By the ergodic theorem 

= -il^^ W"'^) = - /in I /? I d/. < (5.9) 

for 71-almost all t E Vt. For such a t the sequence /3_„(t) exponentially 
converges to zero, and we can consider the following /Z-measurable function 
on Q 

00 

G{t)^-Y.(5^k{t)h{t_k) (5.10) 

k=l 

It follows from (5.8) that for continuous h we have G(t) = U(t), i.e. the 
a priori multi-valued function G turns out to be single-valued. It is the 
necessary condition for solvability of the equation (5.1). 

Remark that this condition is almost sufficient. Namely, if the function 
G{t) is single- valued, then it gives the holomorphic solution in Q* = fl\{0} 
which is continuous up to the boundary. However, this solution has a singu- 
larity at zero since /3(0) = 0. 
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6 Asymptotics of T" along almost all geodesies 
(conditional result) 

In this section we assume that the function U{i) given by (5.5) is not con- 
tinuous up to the boundary. Under this assumption we will calculate its 
asymptotics along almost all (in the sense of harmonic measure) geodesies 
Bq. It turns out to be the following: 

-ln|^(.e-^)|^X^_ (6.1) 
r-i -ln(l-r) Xm 

Remcirk. By Makarov's theorem [31], InZ(i) ~ ln(l — r) for /x-almost all 

geodesies . So, (6.1) coincides with the required asymptotics (1.5) (take in 
account formulas (4.4) and (5.4) for characteristic exponents) 

First let us give heuristic argument yielding (6.1). By (5. 8), (5. 10) we have 

f3^nmit-n) ^ Uit) - Git) (6.2) 
for almost all t = {t,t_i, ...) e H. By (5.9), 



b' 



where the sign ~ means 



lim In ^ = 0. 



Hence, 



Uit.,,) ^ (^)^ (6.3) 



'2 

Setting t^n — r^n^^^-'^^t — re'^, we get r_„ = r''"" ,and hence 

r 1 

6" = In 



(since r is fixed). Comparing this with (6.3), we find 

' — n 
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where «c = Xn{P)/Xfi is the critical exponent written above. 

The main shortcoming of this calculation is that the points t-n don't lie 
near a single radius Bq but are wandering along However, they quite 
often approach almost every geodesic which allows to turn our argument into 
rigorous. 

Let us consider the multi-valued function G{t) on the covering space Q 
given by the series (5.10). This function is correctly defined due to (5.2), and 
satisfies the following equation 



the lift of equation (5.1) on Q. In order to substantiate the asymptotics we 
need some analytical properties of the G. 

Let us put on f2 ~ [0, 1] x S;, measure A, the product of the Lebesgue 
measure on [0,1] and the Bernoulli measure /I on E^. Consider the space 
A(Q,) of functions measurable on D, , analytic on almost all Riemann surfaces 
W{e-) and continuous up to the boundary in their inner topology (see §4). 
Remark that they probably are not defined at 0. 

Lemma 6.1 The function G{t) given by (5.10) belongs to the space A{Q) 

Proof. By the ergodic theorem, /9_„(| t \) ^ e~°'" {n oo) for /I-almost 
all t = (1,^) e r. Applying (4.6) to the function p = — In | /3 |, we get the 
similar asymptotics for all T eW{t) = W{e^). Moreover, for all ^ > there 
exists Q{5, 9) such that 



provided | r |> 5 > 0. Consequentlythe series (5.10) converges uniformly on 
Ws{t) which yields the required statement. g 

There is the Bernoulli measure fx on the solenoids — S5, and one can 
consider the corresponding spaces L'^{Vr)-i k > 0. Now let us prove the 
Main technical lemma. 

Lemma 6.2 For sufficiently small k > and any r e (0, 1], the function 
G{t) belongs to the space L'^(Tr): 



mv{ft)-v{t) 



h{t), 



(6.4) 



P-n{T) \< e^(^) I p_n{t) \< g(^,^)e-(«-^)^ 
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Moreover, for any 5 > the same is true for the function 

Gs{e) = sup I G{r,e) \ . 

r:S<r<l 

Proof. Applying (4.6) to the function p — \n\P\,we get for r > 5: 

\P-n{r,e) |<g(5) I /?_„(!,£) | . 



Consequently, 



Gs{e)\<{\\^\\ooQ{S)rY.\(^-nihe)r, 

n=l 



and further 



f_\Gsrd]i< (II ^^J Q{s)r f: [jp-u r dp- 



But 



J_\ (5.n{t) r dji{t) = j_\ Pnif "t) r dm = J-\ r dm = 1 1 Pn{t) r dm 

Thus, it is enough to show that for sufficiently small k, > the last 
integral exponentially tends to zero as n ^ oo. To this end we will apply 
thermodynamical formalism. 

Let be the Gibbs measure corresponding to the potential —kIh \ j3 |,and 

= Pf{—K,hi I /5 I) be the corresponding pressure. By (4.1), 

I Pn{teo...er,-i) T"- i^K[£o--£n-l] exp(nP«) 

for any point teo - en-i from the cylinder [^o-.-^n-i]- 
As n[£o...en-i] = then 

/ I l'" d/i X i^«[£o--£n-i] exp(nP«) = exp[n{P^ - In 6)]. 

l<Si<n 

Hence, it is sufficient to check that for small enough k, > 

P«-ln6<0 (6.5) 
21 




Figure 4 



Hence, (6.5) holds for any k, G (0, Kq) (Fig. 4), and the lemma is proved. 



So, we have two solutions of the equation (5.1), U{t) = U{t) given by 
(5.4) and G(i) given by (5.10). The first of them comes from the single- 
valued function on VLq, the second is multi- valued on VLq but possesses better 
boundary properties. The difference V(t) = U{t) — G(t) satisfies the homo- 
geneous equation 

mV{Tt)^V{t). (6.6) 

Since U{t) is not continuous up to the boundary (see the beginning of the 
section), while G e A(Q), we have V We are going to find the asymp- 

totics of V along /I— almost all "geodesies" B- = {(r, e) G : < r < 1} as 
r ^ 1. Since G is continuous on almost all sheets, the same asymptotics will 
be valid for U . 
Set 

Vnie) = max I Vir.e) I . 

l/2<r<6/2 ' ^ ' ^ I 

By Lemma 6.2, Vq G L'^. Hence 
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/ In I Vo(s) I djile) < +00. 
Now we can apply the ergodic theorem. It follows 



lim -\n\Vo(a''e) \<0 (6.7) 

for almost all e E T^h. Let us consider a geodesic B-^ corresponding to such 
an e. Let r > 1/2 and n — n{r) is such a number that 

l<r'^< (6.8) 
It follows from (6.6) and (6.7) that 

I V{r,e) H /?„(r,£+) || Vir'^a'^e) \^\ /3„(r,£+) | 

where the sign "a ^ /3" means limMn(/3„/Q;„) < 0. Since by the ergodic 
theorem, | /3„(r, £+) |~ e"" for almost all £+, we conclude 



Prom (6.8) we find 



\im^\nV{r,e)<xM- 



Two last estimates yield 

K < (6.9) 

r^i -ln(l - r) ~ ln6 ^ ^ 

In order to get the opposite estimate let us consider the set 

Xs^{seJ:,:V{l/2,e)>S} 

. As V ^ 0, X{Xs) > for sufficiently small 6 > 0. Hence, almost all 
orbits {cr"e}J^o P^^^ through Xs infinitely many times. For such an e let us 
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consider a sequence n{k) — >■ oo for which a'^^'^h G Xs- Set r„ = (1/2)'' 
Then (6.6) imphes 

I V{rn(k),£) hi Pn{k)irn{k),£+) \ ' \ V {1 / 2 , a^'^'^h) \> 5 \ f3n{k){rn{k),£+) I ■ 

This imphes the inequahty opposite to (6.9). Thus, for almost aU £ G we 
have 

IE '"Hr.^) 

r^i — ln(l — r) m b 

and the required asymptotics (6.1) is proved. 

7 J^" is not continuous up to the boundary 

Let us prove first that one of the derivatives is discontinuous on the interval 
[0, tc] at the real critical point tc (it is a critical point from the thermodynam- 
ical viewpoint; from the dynamical viewpoint, tc is a repelling fixed point ). 
Denote by A = f"'{tc) > 1 the multiplier of this point. 

Lemma 7.1 Let I be a natural number for which A' > 2b. Then /^'^ is 
discontinuous on the interval [0,tc] at tc- 

Proof. Let us hnearize f at tc- 

where A is an analytic function in an neighbourhood of tc, V'(^c) = 0, ij^'{tc) = 
1 (the Konigs function). Clearly, ip can be analytically continued on the 
interval {0,tc], and one-to-one maps it onto the axis (— oo,0]. 
Set 

Then F satisfies the following functional equation: 

^ F {\z)- F {z) ^9 (z). 

Due to the linearization we immediately get the functional equation for 
F : 
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-^F {Xz)-F {z)=9 {z). 

If F is continuous on the semi-axis (—00, 0], then it can be given on it 
by the following series (compare (5.10)): 

which gives the analytical continuation of F through 0. Hence F can be 
analytically continued through tc contradicting Lemma 3.2. g 

Remcirk. The same argument can be applied to any fixed point a e dQ,: 
some derivative of F | Jl° must be discontinuous at a. 

Lemma 7.2 The function F" is not continuous up to the boundary offl. 

Proof.Assuming the reverse, we will show that all derivatives F^") should 
be continuous up to the boundary, contradicting Lemma 7.1. 

Let us consider the functional equations for the derivatives of F (cf. (3.6) 
and (5.3)): 

where hn can be expressed via the derivatives of F of order < n: 



It is convenient to use the following metric in OP: 

d(^,C) = inf/(7) 

where inf is taken over all rectifiable paths 7 connecting z and C,. 

Assume by induction that all derivatives F'^^\ k = 1,2, ...,n, are contin- 
uous in f2 (n=2 is the base of induction). Then the derivatives F^^\k = 
1,2, ...,n — 1, should be Lipschitz continuous with respect to the metric d. 
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By (7.2), hn should possess the same property. But if F^"^^ is continuous in 
fl, n>2, then it can be given by the series (5.10): 

F^^\t) = - V —i-^—^K-iit-k) (7.3) 

for a //-typical inverse orbit t — {to,i-i, •■■}■ It follows from here (taking in 
account the expanding property of / | F) that F^'^'^ is also Lipschitz contin- 
uous with respect to d. Hence, F'^"^^) is bounded in It is enough for 
pin+i) given by the series (5.10) and, hence, to be continuous up to the 
boundary. g 



8 Concluding remarks 

Note that the complex critical exponent calculated in this paper differs from 
the usual real critical exponent a at ^c- In fact, there is the general scheme 
including both cases. Namely, one can associate to any invariant measure v 
on BQq its own critical exponent ai,, i.e. the exponent of power growth of 
an appropriate derivative of the free energy along i/-typical geodesies. More 
specifically, let 

.ln26^ Jn26, 
= [ J -M, a^ = l- { j 

where [a] and {a} denote the entire and the fractional part of a respectively, 
and Xu is the characteristic exponent of u. Then the following general formula 
should be true: 

In I T^"'''Ht) I 
lim ' , „ V = (8-1) 
t-*T,teBr — In l{t) 

for i/-almost all geodesies. 

Indeed, following the scheme of the present paper, we should find the first 
m for which the function (3 = (/')"*/26 has positive characteristic exponent: 

i.e. m > \n2b/xu (assume that l'n.2b/xu is non-integer). Then by formula 
(6.1) we get (8.1). 



26 



For the 5-measure concentrated at the critical point tc we obtain the usual 
formula of the renorm-group theory (see [29] , [13]): 



In 2b 



ciln2 



}• 



ln/'(tc) 



ln/'(tc) 



The similar formula holds for any periodic point on dD, (without any changes 
in the proof). 

The main technical problem in proving (8.1) for general u is related to 
the fact that if m > 2 then the right-hand side of the equation (7.1) for F^'^^ 
is discontinuous on the boundary. The same problem arises if one wants to 
calculate the critical exponent for the measure of maximal entropy on the 
boundary of of the "high-tempcraturc" basin Qi . The true formula should 
be a = 1 — (In 6/ In 2}. It is interesting also to find the complex critical 
exponent ior b — d — 2. 

The other problem is to study the global properties of free energy on the 
whole Riemann sphere. They have to do with Gibbs measures on J(/). 

We finish the paper with the following important remark. One can show 
that the free energy F can be represented as the logarithmic potential of the 
measure of maximal entropy of /. This gives another approach to the circle 
of problems under consideration and a nice relation of the critical exponent 
to the local dimension of the measure of maximal entropy. We are grateful 
to P. Moussa and A. Eremenko for interesting discussions of this point. 
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